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Abstract 

A general continuum-RPA approach is developed to describe electromagnetic transi- 
tions between giant resonances. Using a diagrammatic representation for the three- 
point Green's function, an expression for the transition amplitude is derived which 
allows one to incorporate effects of mixing of single and double giant resonances 
as well as to take the entire basis of particle-hole states into consideration. The 
radiative widths for El transition between the charge-exchange spin-dipole giant 
resonance and Gamow- Teller states are calculated for 90 Nb and 208 Bi. The impor- 
tance of the mixing is stressed. 
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1 Introduction 



Since the very beginning of nuclear physics, electromagnetic transitions be- 
tween nuclear states have been known to be a sensitive tool to investigate 
structure of the nuclear wave functions. In recent years significant theoretical 
efforts have been devoted to describe intensities of the radiative transitions 
between charge-exchange giant resonances (GR) in nuclei [l]-[4]. The interest 
in the problem is related to the possibility to obtain corresponding experimen- 
tal data from the ( 3 He, ^-reaction cross sections. A confrontation of the data 
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with the calculation results would be especially interesting in view of the the- 
oretical speculations on the possible enhancement of the transition intensities 
[l]-[3]- 



It should be mentioned that conclusions on the enhancement have been based 
upon the consideration of the relevant sum rules. However, the main contri- 
bution to the sum rules seems to come from the nuclear states of 2 particle - 2 
hole type (double GR). For instance, the dipole sum rule for the Gamow- Teller 
resonance (GTR) should be mainly exhausted by the giant dipole resonance 
(GDR) built on the top of the GTR, i.e. a double GR, but not by the spin- 
dipole resonance (SDR), a single GR. In this connection the problem of a 
correct description for the mixing of single and double GR arises, because 
it might change significantly the transition intensities. The relevance of this 
problem is indicated by the observation that the calculations of the radia- 
tive transition intensities within TDA [l]-[3] and RPA [4] when taking only 
particle-hole structure of GRs into consideration have strongly underestimated 
the corresponding sum rules. It is noteworthy, that in recent works [5] the im- 
portance of taking into account the GDR admixture has been demonstrated 
to describe correctly the decay amplitude of the first 1~ (two-phonon type) 
state to the ground state. 

In the present work a general expression for the amplitude of the radiative 
transitions between different GRs is obtained for closed-shell nuclei within 
the continuum-RPA using a diagrammatic representation for the three-point 
Green's function. The expression allows us to take the entire basis of particle- 
hole states into consideration as well as to incorporate the important effects 
of mixing of single and double GR. 

We apply the general approach developed in this paper to calculate the inten- 
sity of El transition between the SDR and GT states. The calculation results 
for 90 Nb and 208 Bi nuclei differ markedly from the previous ones [2]- [4], mainly 
due to the effects of mixing of the single and double GR. The corresponding 
experimental data for the latter nucleus is to appear soon [6]. 



2 The amplitude of electromagnetic transitions between giant res- 



We start with the well-known formula for the partial width corresponding to 
a radiative transition between initial |1) and final |2) nuclear states with the 
excitation energies E\ and E 2 , respectively: 
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Here, = Y. a Qk^( x a) is a single-particle multipole operator, (2||Q fe ||l) 
is the reduced matrix element, x is the nucleon coordinate (including spin 
and isospin variables), m is the angular momentum projection, C 7 is an ap- 
propriate dimensional coefficient. In the case of El transitions, considered in 
the present paper to apply the general approach (see Sect. 3), Q i/X = = 

— \er Yxpin) is the isovector electric dipole operator and C 1 = ^ (j^) 
with E 7 — Ei — E 2 being the gamma-quantum energy. 

Let |1) and |2) be states of the particle- hole (p-h) type so that their structure 
can be described within the RPA. Expression for the transition amplitude 
(2||<2fc||l) can be obtained in Green function technique used in the theory 
of finite Fermi-systems [7] to describe the Fermi-system response to a single- 
particle probing operator. Let us consider the diagrams for the 3-point Green 
function (3-point vertex function) describing the transition amplitude of the 
system under the action of the external field Q (Fig. 1). 




Fig. 1. Three-point Green's function describing the amplitude of the transition 
between the RPA states 1 and 2. The solid line presents the nucleon, the dashed one 
the single-particle external field Qk(x), and the wavy lines the transition densities. 

The p-h structure of the state 1 (2) enters these graphs by means of the 
transition fields V\(2){ x ) — I F{ xx ')pi(2){ x ') dx' with pu2)(x) and F(xx') being 
the transition densities and the particle-hole interaction (depicted by the full 
diamonds), respectively. The transition fields are proportional to the residues 

in the poles of the effective operator : Vu2)(w) — — IiHL_MHl ; with V satisfying 

iO — £1(2) 

the equation [7]: 

V(x, u) = V{x) + [ F(xx')A(x'x", uj)V(x", u) dx'dx", (2) 



Here, A(xx',uj) is the free particle-hole propagator, carrying the quantum 
numbers of the appropriate p-h state, the probing single-particle operator 
Vi(2)(a;) has the matrix element M 1(2) = (1(2)|V 1(2) |0) (V = E a V(x a )). 
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The matrix element corresponding to the graphs is: 



(2||Q,||l) = (-) Jl - fe E 

345 \hnk 



x 



n 4 (l - n 5 )(l - n 3 ) 

(e 4 + E 2 - £ 5 ) (£ 4 + E 1 - £ 3 ) (£ 4 + E 2 - £ 5 ) (£ 4 + E l - £ 3 ) 

(1 - n 4 )n 5 n 3 



-<4|h||3>(3||Q fc ||5><5|h+||4> 



_)^+*(4||t;+||5>(5||Q fc ||3>(3|K||4> 

(1 - n 4 )n 5 n 3 \ 

?i-e 3 )J 

n 4 (l - ra 5 )(l - ra 3 ) 



(3) 



(£ 4 -E 2 - e 5 )(s 4 -E x - £ 3 ) (£ 4 - E 2 - £5)^4 - Ei - e 3 ) 



with and £a being occupation numbers and energies of the single-particle 
states, respectively (A is the set of the single-particle quantum numbers, (A) = 
{/j}). The definition of the reduced single-particle matrix elements (A||Qfc||/i) 
by the Wigner-Eckart theorem is taken in accordance with Ref. [8]. 

The Eq. (3) can be transformed to express the transition amplitude in terms 
of the RPA X and Y amplitudes (see, e.g., Refs. [5,9]): 



(2110*111) = v /(2j 1 +i)(2j 2 +i)5:(-) ji+i+ ^ 3 \ j2k Jl \ 

345 J3J4J5 



X 



(-) Jl+j2+fc (5||Q fc ||3> (iSJlg 1 



\Ji- 



-<3||Q fc ||5> (X$X, 



(2) , 
45 + 



+J2+j3+j5v( 1 )v( 2 )\ 

J- 43 1 45 ) 
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with the use of the following representations for the X and Y: 



21 " ( j V2 Ji +T (iTT E X - e 2 ) {b) 

21 " ( j V2 Ji + 1 (.1 - - (b) 



One sees in Fig. 1 that there is an asymmetry between the vertex 1,2 ("dressed") 
on the one hand and vertex 3 ("undressed") on the other. This is because we 
have not incorporated effects of particle-hole interaction in channel 3, or, in 
other words, we neglected the effect of the virtual excitation of the corre- 
sponding giant multipole resonance. This effect has been well-known to play 
a crucial role in describing the neutron radiative capture since the model of 
direct-semidirect capture was suggested [10] (see also Ref. [11] for more recent 
developments). This additional contribution which takes into account virtual 
excitation of the giant resonance is shown in Fig. 2. 
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Fig. 2. The generalized three-point Green's function which takes into account the 
virtual excitation of multipole giant resonance in the channel 3. 



The wavy line in the channel 3 depictures the total p-h propagator A describing 
off-shell excitation of all p-h states (including the corresponding GR) with 
appropriate spin and parity. Using the fact that AQ = AQ, one can see that 
the sum of the graphs of Figs. 1-2 is equivalent to the graphs of Fig.l only 
if one replaces the external field Qk(x) by the corresponding effective field 
Qk(x,u) satisfying Eq.(2). 

The Eq.(3) with the substitution Qk{x) — > Qk{ x ,uj) can be transformed fur- 
ther with the use of the spectral decomposition for the Green's function of 

the single-particle radial Schrodinger equation gt X ) ( rr ' ? e ) — E — 1 
to include explicitly the single-particle continuum: 



|C fc ||i> = (-) 



Jl-fc 



E ((-) Jl+j2+fc G 4 (E 1 , E 2 ) + <%(-E u -E 2 j) 



E^3^5 ((-) Jl+J * +k G 3 , 5 (E u E 2 ) + Gl 5 (-E l} -E 2 ] 

3,5 



(7) 



G 4 (E U E 2 ) = E *(3)(4)(5)(4|uJ (g(5)(e 4 + E 2 ) -Y! n 5 

|3)(3| 



(3)(5) \ 5 S± + E 2 -e 5/ 



*Qk 5(3) (£4 + ^l) - E ^ll 4 )> 

\ 3 ^4 + -ca — £3 / 

Gsfi^Eh) = E*(3)(4)(5)- f^*' 3) F (3h(g(4)(ea - #i) - <7 W (e6 - ^)K|5), 
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J 2 k Ji\ 



(3)(4)(5) 



L ((4)lk 2 + ||(5))((5)||g fe ||(3))((3)||t; 1 ||(4)), 



where we have decomposed the reduced matrix elements (A||w||/i) = ((A) ||t> || (fi))(X\v \fi) 
into its spin-angular ((A)||v||(/x)) and radial (\\v\fi) = J Xx( r ) v ( r )x^( r ) dr 
parts. The radial integrals (for brevity omitted) can be easily recovered, for 
instance, 



(4|u^ ( 5)(e 4 + E 2 )Q k g (3) {e A + £?i)ui|4) = 
/ XA(r)vi(r)g {5) (rr',e 4 + E 2 )Q k (r')g {3) (r'r",e 4 + E 1 )v 1 (r") Xi (r") drdr'dr". 

The superscript T means transposition and Yl runs over single-particle states 
with the same I and j. 

The Eq. (7) gives the general expression for the amplitude of the electromag- 
netic transitions within the CRPA. Apart from the single-particle continuum, 
it contains also contribution from the double phonon configurations. To see 
this, let us decompose the amplitude into "direct" and "semi-direct" parts ac- 
cording to the partition Qk = Qk + {Qk — Qk)- The first, direct, part describes 
the amplitude in the single-phonon space and it has been only taken into ac- 
count in all previous theoretical considerations [l]-[4]. As for the semi-direct 
part, the following approximate representation holds: 

Qk{x^)-Qk{x)^Y.^^- (8) 



with qi(x) being the transition fields for the doorway states forming the GDR. 
Thus, one can see, that the semi-direct part originates from admixtures of 
double-phonon configurations, like \2) + ® The contribution to the 

transition amplitude is determined by both the matrix elements Mj = (i\Qk\0) 
(states |1) or |2) serve as spectators in these electromagnetic transitions) and 
the admixture amplitudes ctj. 

It can be seen from Eqs. (7,8) that the mixing is described within first order 
perturbation theory with respect to the particle-hole interaction. Neverthe- 
less, the approximation is satisfied in most cases due to the small values of 
the mixing amplitudes, that can be also seen from the calculated transition 
amplitudes being much less than the amplitude between the GDR and the 
ground state. 
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3 El transitions between the spin-dipole and Gamow- Teller reso- 
nances. 



As an application of the approach described in the previous section, we cal- 
culate intensities and branching ratios for El decay of the charge-exchange 
spin-dipole resonance to some Gamow- Teller states (GTR and low-lying satel- 
lites) identified in Ref. [12]. To calculate the GT and SD strength distribu- 
tion, we use the CRPA approach of Refs. [4,13,14]. The phenomenological 
isoscalar part of the nuclear mean field and the zero-range Landau-Migdal 
particle-hole interaction are the ingredients of the approach along with some 
selfconsistency conditions. The approach has allowed one to reproduce very 
well the experimental GTR energies in 90 Nb and 208 Bi nuclei [14]. For each 
spin Jsdr = (0~> 1) °f the SD components we choose just one state with the 
maximal strength to be the SDR, lying at (26.6, 24.1) MeV in 208 Bi and at 
(22.4, 20.1) MeV in 90 Nb, and for J£ DR = 2~ we choose three states with the 
maximal excitation energies, lying at (20.0, 21.3, 23.1) MeV in 208 Bi and at 
(13.5, 14.4, 17.4) MeV in 90 Nb. 

The details of the calculation of the effective fields for the SD and GT cases 
can be found in Refs. [4,13,14], and for the isovector dipole one - in Ref. [11]. 
In the latter case we take the velocity dependent particle-hole forces along 
with the isovector Landau-Migdal interaction that allows us to reproduce in 
the calculations both the experimental GDR energy and the observed excess 
of the energy-weighted sum rule over the TRK one. Because some transition 
energies are rather close to the GDR energy, the influence of the GDR spread- 
ing becomes important. We use the substitution u — > u + \T{u)/2 in Eq. (2) 
to simulate the effect of GDR spreading (the mean energy- dependent doorway 
spreading width T{u) has been chosen as in Ref. [11]). 

We calculate the reduced matrix elements (GTR\\T))~\\SDR, J) according to 
Eq. (7) and Eq. (3) (with and without taking into account the virtual GDR 
excitation, respectively), and then the corresponding radiative widths T 7 and 
T 9 according to Eq. (1). The calculated widths for the El decay of the above- 
mentioned SD states to the GT states at the excitation energy E GT are listed in 
Table 1. The corresponding branching ratios are calculated as the ratio T 7 /r J - 
with being the mean doorway-state spreading width chosen to reproduce 
the experimental total SDR width in the strength-function calculations (1^ = 
4.7 MeV for 208 Bi [13], T l = 4. MeV for 90 Nb). 

As can be seen from the calculation results listed in the Table 1, taking into 
account of the virtual GDR excitation affects the results drastically. In the 
cases where the transition energy is well below the GDR energy, a destructive 
interference between the direct and semidirect parts of the amplitude always 
occurs, in agreement with both the qualitative description in terms of the ef- 
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Table 1 



Calculated radiative widths and branching ratios for the El transitions between 
SDR and GT states in 90 Nb and 208 Bi nuclei without (T°,b°) and with (r 7 ,6 7 ) 
taking into account the virtual GDR excitation, respectively. 





E GT , MeV 




r°, keV 






T 7 , keV 




Nucleus 






(6°, xlO- 


4 ) 


(6 7 , xl0~ 4 


) 




/Jsdr 


O - 


1- 


2 _ 


O - 


1- 


2 - 




8.9 


3.40 


1.03 


0.10 


3.68 


0.02 


0.008 


90 Nb 




(8.5) 


(2-6) 


(0.3) 


(9.2) 


(0.05) 


(0.02) 




0.0 


0.21 


0.10 


0.08 


1.62 


0.13 


0.37 






(0.5) 


(0.2) 


(0.2) 


(4.0) 


(0.3) 


(0.9) 




15.5 


2.28 


0.70 


0.09 


6.49 


0.03 


0.008 






(4.9) 


(1.5) 


(0.2) 


(13.8) 


(0.05) 


(0.02) 




9.3 


0.15 


0.03 


0.14 


0.80 


0.43 


0.22 






(0.3) 


(0.05) 


(0.3) 


(1.7) 


(0.9) 


(0.5) 




8.1 


0.24 


0.02 


0.04 


1.05 


0.26 


0.16 






(0.5) 


(0.03) 


(0.08) 


(2.3) 


(0.5) 


(0.3) 


208 Bi 


6.4 


0.23 


0.007 


0.001 


0.58 


0.04 


0.03 






(0.5) 


(0.01) 


(0.002) 


(1.2) 


(0.1) 


(0.06) 




4.9 


1.05 


0.0002 


0.002 


2.36 


0.003 


0.01 






(2.2) 


(0.0005) 


(0.003) 


(5.0) 


(0.007) 


(0.03) 




4.5 


0.50 


0.0003 


0.001 


1.24 


0.004 


0.005 






(1.0) 


(0.001) 


(0.003) 


(2.6) 


(0.01) 


(0.01) 




3.8 


0.58 


0.004 


0.006 


1.19 


0.02 


0.02 






(1.2) 


(0.01) 


(0.01) 


(2.5) 


(0.03) 


(0.04) 




1.2 


0.1 


0.19 


0.02 


0.21 


0.42 


0.04 






(0.2) 


(0.4) 


(0.04) 


(0.4) 


(0.9) 


(0.1) 



fective charge [8] and previous quantitative consideration of Ref. [5]. In the 
cases when the transition energy is close to the GDR energy or exceeds it, 
the opposite situation takes place and one can see a significant enhancement 
of the corresponding amplitudes. This circumstance along with the steep 
dependence of the El radiative width leads to the fact that the widths for 
SDR decay to some low-lying GT states can be comparable to those for SDR 
— > GTR El transitions. Some of the calculated branching ratios are of order of 
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10~ 3 , that seem to be accessible in experiments. In fact, the statistical back- 
ground hinders detection of SDR — ■> GTR decay, whereas it is small enough to 
allow detection of SDR decay to low-energy GT states which has been already 
observed in ( 3 He,t7) experiments on 208 Pb [6]. For instance, strong gamma 
transitions have been observed to the unresolved GT states with E=3.8, 4.5, 
4.9 MeV [6]. The energy of these gamma transitions (about 14 MeV) is very 
close to the energy of the GDR, and the predicted branching ratio (~ 10~ 3 ) 
can explain the observed transitions. However, we have to mention that the 
present calculations overestimate the mean energy of the SDR by about 2 
MeV [13]. Bearing in mind strong energy-dependent renormalization of the 
transition amplitude found in this paper, this can lead to a noticeable differ- 
ence between the calculated and observed branching ratios. We have not also 
taken the finite widths of the GTR and the SDR into consideration. 

Special interest for testing the model lies in the El transition from the GTR 
to a 2.9 MeV 2~ SD state in 208 Bi, which was observed recently at KVI [6]. 
The latter state has presumably a (vrl/ig^X^l^is^) -1 particle-hole structure 
and has a predicted excitation energy of 2.4 MeV. For this El transition we 
have found a large branching ratio 6 7 = 4.3 ■ 10~ 3 due to the fact that the 
transition energy almost exactly coincides with the GDR energy, therefore we 
obtain maximal enhancement. 

Finally, it is also worth to mention the sensitivity of the radiative widths to 
the spin of the SDR components due to their different energy position that 
could be in principle used to get experimentally information on their energy 
splitting. 



4 Conclusion 



A general continuum-RPA approach has been developed in the present pa- 
per to describe the intensities of electromagnetic transitions between different 
giant resonances. We have used a diagrammatic representation for the three- 
point Green's function to derive the expression for the transition amplitude. 
The expression has allowed us to take the entire basis of particle-hole states 
into consideration as well as to incorporate important effects of mixing of 
single and double giant resonances. The radiative widths for El transition be- 
tween the SDR and GT states (including GTR) have been calculated for 90 Nb 
and 208 Bi nuclei. Some of the widths seem to be large enough to be accessible 
in experiments. 
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